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A nematic liquid-crystal gel is a macroscopically homogeneous elastic medium with the rotational 
symmetry of a nematic liquid crystal. In this paper, we develop a general approach to the study 
of these gels that incorporates all underlying symmetries. After reviewing traditional elasticity 
and clarifying the role of broken rotational symmetries in both the reference space of points in 
the undistorted medium and the target space into which these points are mapped, we explore the 
unusual properties of nematic gels from a number of perspectives. We show how symmetries of 
nematic gels formed via spontaneous symmetry breaking from an isotropic gel enforce soft elastic 
response characterized by the vanishing of a shear modulus and the vanishing of stress up to a 
critical value of strain along certain directions. We also study the phase transition from isotropic to 
nematic gels. In addition to being fully consistent with approaches to nematic gels based on rubber 
elasticity, our description has the important advantages of being independent of a microscopic model, 
of emphasizing and clarifying the role of broken symmetries in determining elastic response, and 
of permitting easy incorporation of spatial variations, thermal fluctuations, and gel heterogeneity, 
thereby allowing a full statistical-mechanical treatment of these novel materials. 
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I. INTRODUCTION 

The term liquid crystal has traditionally 

been used to describe phases of matter that exhibit 
anisotropies characteristic of crystals but that under ap- 
propriate conditions flow like a liquid. These phases typi- 
cally have symmetries intermediate between that of a ho- 
mogeneous isotropic fluid and that of a three-dimensional 
periodic crystalline solid. Indeed one can provide an 
almost complete characterization of a liquid-crystalline 
phase by specifying its symmetry. For example, the ne- 
matic phase, which is spatially homogeneous yet optically 
uniaxial has -Dooh symmetry. The typical phase sequence 
for a thermotropic liquid crystal on cooling begins with 
an isotropic fluid and ends with a crystalline solid after 
passing through nematic, layered smectic-^4 and smectic- 
C, and possibly hexatic phases. 

There is, however, a large variety of materials that 
have the same macroscopic symmetry as fluids, but that 
cannot flow: they are macroscopically homogeneous and 
isotropic elastic media with a nonvanishing shear modu- 
lus that provides resistance to shear distortions. We will 
refer to these materials, which include everything from 
glasses to elastomers or rubbers ||, as gels |Q. One can 
imagine phases arising from a reference state of a gel 
(rather than a liquid) with the same macroscopic sym- 
metries as conventional liquid crystals. As we shall dis- 
cuss more fully below, these phases do in fact exist ||- 
jL2|, and, because they cannot flow, they have mechanical 
properties and mode structures that differ significantly 
from those of standard liquid crystals. We will call these 



phases "liquid-crystal gels" because they are gels with 
the symmetry of conventional (fluid) liquid crystals. In 
this paper, we will develop a powerful and general formal- 
ism to describe nematic gels and use it to explore their 
remarkable properties. We will focus particularly on ne- 
matic gels that form via spontaneous orientational sym- 
metry breaking from an isotropic gel phase. Our formal- 
ism can be generalized to treat other liquid-crystalline 
gel phases. 

There are a large number of experimental realizations 
of liquid-crystal gels. Of particular interest to us are 
liquid-crystal elastomers ||- These materials, which 
are formed by weakly crosslinking either side-chain j^] 
or main-chain j|] polymers, combine the enormous ex- 
tensibility of rubbers with the orientational properties 
of liquid crystals. They are, therefore, of considerable 
technological importance. The existence of the rub- 
bery crosslinked network appears to have relatively lit- 
tle effect on liquid-crystalline phase behavior, and the 
standard thermotropic nematic, cholesteric, smectic-A, 
and smectic-C phases have their elastomeric counterparts 
P,^3[. The elastic properties of these phases do, however, 
crucially depend on whether a given liquid-crystalline or- 
der was established before or after crosslinking. Liquid- 
crystal gels can also be prepared in other ways, for ex- 
ample by polymerization of monomer solutes in a liquid- 
crystalline solvent Jl0[ , or by confining conventional liq- 
uid crystal inside a dilute flexible matrix, such as, e.g., 
aerosil 0Q. 

To fully characterize liquid-crystal gel phases, two 
complementary basic questions must be addressed: (1) 
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What effect does liquid-crystal order have on the gel elas- 
ticity? (2) How does the rigidity of the underlying gel af- 
fect liquid-crystal order and its stability to fluctuations? 
Our primary concern in this paper will be with question 
(1) applied to nematic gels. With regard to question (2), 
gel elasticity of weakly crosslinked elastomers has rela- 
tively little effect on the existence of liquid-crystalline 
phases. On the other hand, liquid-crystal elastomers ap- 
pear experimentally to be more strongly ordered than 
their fluid counterparts. For example unlike conventional 
nematics, which are milky and therefore have domain 
sizes at the scale of visible light, elastomer nematics are 
clear, indicating orientational order extending beyond a 
micron. This property directly and clearly follows from 
our model of nematic elastomers, as well as early work 
jwj], and it will be explored in more detail in a future 
publication. [ 51| 

The strong interplay between broken symmetry and 
the nature of long-wavelength excitations of ordered 
phases is a major theme of physics . Symmetry prin- 
ciples dictate that ordered thermodynamic phases that 
break a continuous symmetry have low-energy distor- 
tions, or "soft" modes, that are described by an elastic 
energy depending only on gradients of these Goldstone 
fields whose spatially uniform increments take the system 
to symmetry-equivalent states. The form of this elastic 
free energy is uniquely determined by the properties of 
the reference phase whose symmetry is broken and by 
the nature of the broken symmetry itself. Conventional 
nematic liquid crystals break the rotational isotropy of 
an isotropic homogeneous liquid, and they are charac- 
terized by the Frank elastic free energy which is a 
functional of the Goldstone field n, the Frank director 
specifying the direction of molecular alignment. As il- 
lustrated in Fig. |l], nematic phases of liquid-crystal gels 
that form from an isotropic gel state also spontaneously 
break rotational isotropy. Their long-wavelength elastic 
energy, however, differs significantly from the Frank free 
energy of conventional nematics because the reference gel 
state (unlike a reference fluid state) has a nonvanishing 
shear modulus. The elastic energy of a spontaneously 
formed nematic gel was first calculated by Golubovic and 
Lubensky (GL) [0 in their study of a model isotropic 
elastic medium that undergoes a phase transition to a 
uniaxial state when its shear modulus becomes smaller 
than a critical value. They found that the Goldstone 
fields of a uniaxial gel are displacement fields and that 
their associated elastic energy is expressed as a function 
of standard strains of a solid. Normally, the elastic en- 
ergy of a uniaxial elastic medium is characterized by 5 
independent elastic constants. A nematic gel that forms 
spontaneously from an isotropic gel is significantly softer 
than a conventional uniaxial solid: it is characterized by 
"soft" elasticity in which the elastic constant C5 associ- 
ated with shears in the plane containing the anisotropy 
axis vanishes and in which stress vanishes up to critical 
values of certain shears [Q- piifl . 

Many of the properties of nematic elastomers can be 



explained by an elegant and remarkably simple exten- 
sion of standard rubber elasticity [£2| in which ne- 
matic order leads to an anisotropic step-length tensor 
for random-flight polymer segments between crosslinking 
points. This "neoclassical" theory of rubber elasticity, 
which describes a particular realization of an anisotropic 
gel exhibits "soft" elasticity [^3|,[l4) in accord with the 
general symmetry-based predictions of GL. 
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FIG. 1. A cartoon of a liquid-crystal gel undergoing an 
Isotropic-Nematic transition, accompanied by a spontaneous 
uniaxial distortion. 

In this largely pedagogical paper, we will explore the 
elastic and orientational properties of nematic gels from 
a perspective that is an extension of that of Ref. |l7]] 
rather that of rubber elasticity JtJ . In particular we will 
describe elastic properties mostly in terms of nonlinear 
strain tensors familiar from the elastic theory of solids 
and membranes |^J,^5| rather than the perhaps more 
fundamental Cauchy deformation tensor used in rubber 
elasticity ]23] from which the nonlinear strain tensors can 
be constructed. We will usually measure strain relative to 
some equilibrium state in which the strain is zero rather 
than relative to a state at the time of preparation as is 
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common practice in the theory of (incompressible) rub- 
bers. Our approach based on nonlinear strains that are 
invariant under arbitrary rotations of the sample (or, al- 
ternatively, as we shall see, the reference state) allows us 
to keep track of rotational invariances with relative ease. 
It is particularly well suited, as we will show in a future 
publication p6| , to the treatment of renormalized elas- 
ticity arising from the interplay of thermal fluctuations 
and nonlinear elasticity. The approach is also convenient 
for the discussion of external-field-induced instabilities of 
an equilibrium phase. Most importantly, our formalism 
elucidates the origin of the novel soft elasticity of nematic 
gels, making it clear that it arises from general symme- 
try principles common to any spontaneously uniaxially 
ordered elastic medium and is not limited to any specific 
model of such materials. Our description has the disad- 
vantage compared to the rubber elasticity approach that 
it does not naturally treat the very large (as much as 
400% j27|| ) extensions that can arise in elastomers. 

Though gels are macroscopically isotropic and homo- 
geneous, they are always randomly anisotropic and in- 
homogeneous at sufficiently short length scales. Conse- 
quently, there is a local preferred direction of orienta- 
tional (and spatial) order that acts as a random orienting 
(and pinning) field. These quenched fields are certainly 
present in anisotropic gels: optical observations in thin 
films provide direct evidence of their existence p8| , p9[ . 
One consequence of such random quenched local fields 
is that elastomers crosslinked in the isotropic phase and 
cooled into the nematic phase exhibit a polydomain ori- 
entational structure, which disappears at a polydomain- 
monodomain transition when a sufficiently large external 
stress is applied 30 1. Such random static fields can eas- 
ily be incorporated in our formulation of nematic gels. 
Their study is in principle necessary to understand com- 
pletely the effect of gel matrix on nematic order [question 
(2) above]. Such an investigation would parallel a body 
of work on conventional liquid crystals confined inside 
the quenched, random, but (nearly) nondeformable envi- 
ronment of rigid gels, such as e.g., an aerogel. ]TT| , ^l| , ^2[ 
Experience with these rigid systems indicate that ran- 
dom fields might become qualitatively important at suf- 
ficiently long scales. Nevertheless, in this article, we will 
ignore completely effects of random fields and concen- 
trate on properties of anisotropic gels formed from ideal 
isotropic homogeneous gels. By focusing on gels in which 
crosslinks are dense and well-percolated, we will also have 
nothing to say about the nature of the vulcanization tran- 
sition itself fl33fl . 

We will also leave for the future |34| the analy- 
sis of "semi-soft" nematic elastomers ]I8|J20| , [35| that 
are prepared by polymer crosslinking in the nematic 
phase. These materials are characterized by a small non- 
vanishing elastic modulus C5 and nonlinear stress-strain 
curves with a small but nonvanishing stress up to large 
strains. 

This paper is organized as follows. Section II reviews 
the standard Lagrangian theory of elasticity and estab- 



lishes notation for sections that follow. It introduces the 
reference space consisting of points in the undistorted 
medium and the target space into which these points 
are mapped. This section emphasizes the two distinct 
rotational invariances of isotropic elastic media, namely 
invariance with respect to rotation of the deformed sam- 
ple itself (rotations in the target space) and invariance 
with respect to rotation of points in the original refer- 
ence material that map to particular points in the target 
space (rotations in the reference space). Section II also 
discusses the standard nonlinear strain tensor, the right 
Cauchy-Green strain tensor, that is invariant with re- 
spect to rotations in the target space and introduces an 
alternative nonlinear strain tensor, the left Cauchy-Green 
strain tensor, that is invariant with respect to rotations 
of the reference space but transforms like rank-2 tensor in 
the target space. Section III elaborates on the model con- 
sidered in Ref. (Tt]]. It shows in particular that the elas- 
tic energy of the anisotropic gel phase expanded only to 
harmonic order in the nonlinear strain does not preserve 
the rotational invariance of the original energy with re- 
spect to rotations in the reference space. It also discusses 
the isotropic-to-anisotropic transition in terms of the al- 
ternative strain tensor that preserves rotational invari- 
ance in the reference space. Section IV discusses a model 
with both strain and the symmetric-traceless tensor or- 
der parameter Qij of the nematic state. It shows that the 
"soft" elasticity of the anisotropic (nematic) state arises 
because the nematic order parameter can relax strains in 
the plane containing the anisotropy axis, as first shown 
by Olmsted Section IV also derives the elastic en- 
ergy for nematic glasses deep in the ordered phase, where 
biaxial fluctuations can be neglected and nematic prop- 
erties can be described completely by the Frank director. 
This theory is expressed in terms of generalized nonlin- 
ear strains that are functions of strain and relative di- 
rector orientation and that are invariant with respect to 
arbitrary simultaneous rotations of the director and mass 
points. Section V explores the relation between the the- 
ory presented here and the neoclassical elastomer theory. 
Section VI concludes with a discussion of some of the 
many interesting open problems, such as instabilities of 
elastomers induced by external perturbations (e.g., elec- 
tric or magnetic fields) and the effects of thermal fluctu- 
ations and quenched local random anisotropy fields, that 
can conveniently be addressed through our formulation 



II. CLASSICAL LAGRANGIAN ELASTICITY 

Classical elasticity |24j provides a phenomenological 
description of the energy associated with slowly |36f vary- 
ing distortions of an elastic body from its equilibrium 
configuration. As discussed in the Introduction, it is a 
symmetry-restricted theory of the low-energy Goldstone 
modes associated with spontaneous translational symme- 
try breaking. In this section, we will review the classical 
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theory of Lagrangian elasticity |37| , introducing concepts 
that will be important for our study of spontaneously 
uniaxial nematic elastomers. 



A. Strain 

The equilibrium unstretched medium occupies a re- 
gion of a Euclidean 3-space, which we will call the refer- 
ence space Sr. Mass points in this medium are indexed 
by their vector positions x = (xi, X2, X3) = {x,y,z) in 
Sr, which are their positions in the unstretched medium. 
When the medium is distorted, the point originally at x 
is mapped to a new point R(x) = (i?i(x), i?2(x), i?3(x)) 
in Euclidean space. We will refer to the space of points 
defined by R as the target space St- Since there is no 
distortion when R(x) = x, it is useful to introduce the 
displacement vector u(x) that measures the deviation of 
R from x: 

R(x) = x + u(x). (2.1) 

Both Sr and St are Euclidean, with distances deter- 
mined by the unit metric: dx 2 — dxidxi, and dR 2 — 
dRidRi, where the Einstein summation convention on 
repeated indices is understood. j3£| It is often interesting 
to consider generalizations of the above picture to a D- 
dimensional reference space and a d > D target space, for 
example to describe D = 2-dimensional tethered mem- 
branes fluctuating in d = 3-dimensional real space [ p5[ . 
In this paper, however, we will restrict our attention to 
D = d = 3, leaving discussion of membranes to a future 
publication. [p6| 

Distortions that vary slowly on a scale set by micro- 
scopic lengths of the reference material (interparticle sep- 
aration in a glass, distance between crosslinks in an elas- 
tomer, etc.) are described by the Cauchy deformation 
tensor ||@, 

OR 

A« = ^ = d J R l = % + na , (2.2) 

where 

rjij = djUi (2.3) 

is the displacement gradient tensor. Throughout the pa- 
per, we will often use matrix notation in which M is the 

matrix with components My and M T is the transpose 
matrix with components Mj%. 

The energy of the distorted state relative to the undis- 
torted one depends on the how much the target space 
is stretched relative to the reference space, i.e., by how 
much the distance between two nearby points changes in 
the mapping from the reference to the target: 

dR 2 — dx 2 = 2uijdx{dxj , (2-4) 

where 



= \ (AfcjAfcj - Sij) or g = I (A T A - £) 
= j{d l u j + djUi + diUkdjUk) 

= §(% + Vji + VkiVkj)- (2.5) 

is the familiar nonlinear Lagrangian strain tensor pi[ , 
also called the right Cauchy-Green strain tensor or simply 
the Green strain tensor |39| . It is symmetric by construc- 
tion. It is also invariant (lie., transforms as a scalar) un- 
der arbitrary rotations of the target space vector R, i.e., 
if Ri is replaced by R^ = OTij Rj , where Orij is an arbi- 
trary rotation matrix, iiy does not change. On the other 
hand, Ujj transforms like a rank-2 tensor under rotations 
of the reference space, i.e., if cc, — » x[ = O^Xj, then 

l^yg 1 . (2.6) 

Isotropic solids, e.g., the glasses and gels of interest to 
us, are (statistically) invariant under arbitrary rotation 
O r in the reference space Sr. Crystals, on the other 
hand, have lower symmetry and are invariant only under 
a point subgroup of all rotations O . 

In contrast, invariance with respect to arbitrary rota- 
tions in St is a property of all elastic media in the absence 
of external aligning fields, whether they be isotropic, 
crystalline, or wildly inhomogeneous. Thus, because it 
by construction incorporates the Ot invariance, in most 
instances, « y - is the strain tensor in terms of which elastic 
theory is most conveniently formulated. However, here 
we are interested in systems (gels) that exhibit rotational 
invariance in the reference space, i.e., an Or invariance 
of Sr, and, therefore, a distinct left Cauchy-Green strain 
tensor, 

Vij = \ (Aik^jk - S^) or g = I (AA T - £) 
= \{diUj + djUi + d k Uid k Uj) 
= |(%' + Vji + VikVjk), (2.7) 

is useful p| . This tensor is invariant under arbitrary 
rotations O r in Sr, but it transforms like a rank-2 tensor 
under rotations 0_ T in St- 

In what follows, we will simply refer to u and v as right 
and left strain tensors, respectively. 

The left strain tensor v can be contracted with other 
target-space tensors, such as the Maier-Saupe-de Gennes 
nematic order parameter Qij, or the electric field Ei, to 
form scalar invariants such as TrvQ or EiVijEj. In con- 
trast, the contractions TruQ and EiUijEj are not scalars 

since u does not transform like a tensor in the same space 
as the tensors Qij and EiEj. In the absence of external 
aligning fields such as E that effectively render the tar- 
get space anisotropic, the right strain u provides a com- 
plete description of elastic distortions, even if, as is the 
case for crystals, the reference space is anisotropic. If 
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Sr is isotropic and there are external fields breaking the 
isotropy of St, then u cannot provide a similar complete 
description, but the left strain v can. On the other hand, 
the left strain cannot provide a complete description if 
the reference space is not isotropic. For example, semi- 
soft elastomers crosslinked in the nematic phase with a 
director no, which specifies a direction in Sr, are in- 
variant under the simultaneous rotations of no and x, 
no — > O^Uq and x — > 0~ x, but not under rotations, 

x — > QZ^x, of x alone. The left strain v is a scalar in 
Sr, and it cannot be contracted with the reference space 
vector no. Thus, it is impossible to construct scalar in- 
variants involving v and no and to construct a free energy 
in terms of v that reflects the anisotropy of Sr. If Sr 
is anisotropic and there are external fields breaking the 
rotational invariance of St, then only the deformation 
tensor A can provide a complete description of the en- 
ergy of elastic distortions. 

B. Isotropic systems 

For most gels, the reference space is macroscopically 
isotropic and homogeneous, i.e., like an isotropic fluid, it 
is invariant under x — > T + 0~ X x for arbitrary transla- 
tions T and rotations O in Sr. Thus, the elastic energy 

is invariant under R(x) — * O r R(T + 0~ x). The invari- 
ance under rotations 0_ T of the target space is easy to 
understand: different physical orientations of the mate- 
rial (even if arbitrarily distorted) have the same energy. 
Invariance under O r is somewhat more subtle though 
complementary. Figure provides a useful graphic repre- 
sentation of this invariance in two dimensions. Consider 
a circle of radius r in the reference space consisting of 
the points x = r(cos fa sin <ft) = (r,fa. Under distortion, 
it is mapped onto some closed curve in St consisting of 
points R(</>). Thus, the point (r, fa) in Sr is mapped 
to the point Ri = R(0i) in St, (f, fa) is mapped to 
R2 = H(fa), and so on. Under a rotation through 9 in 
Sr, <j> — > + 9. Because of the isotropy of the undis- 
torted, reference state, the energy is not changed if the 
points (r, <p — 9) rather than (r, fa are mapped to the 
points R(</>), i.e., if (r, fa — 9) is mapped to Ri, (r, fa — 9) 
to R2 and so on. 

Care must be taken to incorporate the above symme- 
tries in the free energy density of such homogeneous and 
isotropic gels. Invariance with respect to translations T 
in Sr is enforced by requiring that the free energy density 
depend only on spatial derivatives of R with respect to 
x, i.e., depend only on A and possibly higher derivatives 
of R. Under rotations in Sr and St, the Cauchy strain 
tensor Ay transforms according to 

A y ^O m ||g=0 Tiib A«0^.. (2.9) 

The free energy density / of an isotropic gel is invariant 
under independent Q_ T and O r rotations and must satisfy 



/(A) = /(a T A£ i? 1 )- (2.10) 
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FIG. 2. Schematic representation of mappings from the 
reference space Sr to the target space St- The points 
xi = (r, <f>i) and X2 = (r, fa) in Sr are mapped, respectively, 
to the points Ri and R2 in St- There is a strain energy E$ 
associated with this mapping. For isotropic reference spaces, 
mapping of points in Sr, first rotated by 6 inside the refer- 
ence space (described by a rotation matrix O^ 1 ) to the same 
set of points in St, i.e., mapping points x' x = (r, <f>i — 8) 
and x 2 = (r, fa — 0) to Ri and R2 clearly produces the 
same energy Es as the unrotated mapping. Subsequent tar- 
get-space rotation O t by Qt of the resulting distorted state 
with points Ri and R2 mapped to R' x and R 2 costs no energy. 
The transformation R(x) — » Q_ t R(CX r 1 x) between energeti- 
cally-equivalent nematic states is the Goldstone mode respon- 
sible for the novel elastic properties of nematic elastomers. 

It must, therefore, be constructed from the scalar in- 
variants Tr(AA T )" and detAA_ T = (det A) 2 @. Alter- 
natively, the free energy can be equivalently expressed in 
terms of u or v, with the respective invariances 



/y=/(Q Ri a- 1 ) 

f(g) = f(Q r gQ^ 1 ), 



(2.11a) 
(2.11b) 



which are enforced by allowing only fully contracted pow- 
ers of strain tensors to appear. The energies f(u) and 
f(v) can be derived from /(A) using 

det AA_ T = exp Tr ln(| + 2^) 
= expTrln(|+ 2^), 

T m n = Tr [i (A T A ~£)] n = Tr [\ (AA_ T - |)] " 
= Trg". (2.12) 

Thus, f(u) and f(v) depend only on Tru™ = Tru™, and 
f(v) is the same function of v that f(u) is of u. For 
the discussion of encoding these two Or and Ot symme- 
tries at the harmonic level in the phonon variable u, see 
Appendix ^|. 

Although many of the properties of nematic elastomers 
follow directly from the above invariances, in what fol- 
lows, it will be useful to have explicit forms for the elastic 
free energy density. A model free energy in terms non- 
linear strain tensor u, up to fourth order in u is 



f(u) = ±A(Tr^) 2 + ^Tr^ 2 - CTr^ 3 
+L>'(Tru 2 ) 2 - E'TruTru 2 . 



(2.13) 



As just discussed, this free energy can equally well be 
expressed in terms of v merely by replacing u by v. In- 
variances with respect to rotations in Sr and St are en- 
forced in f(u) in different ways. Symmetry under O t is 
enfo rced by the construction of the strain tensor u, Eq. 
(2.5), which, being a scalar in St, is automatically in- 



variant under O m . Invariance under O „ is enforced by 



only allowing terms in f(u) that transform as a scalar 
under O , i.e., only fully contracted powers of u. In con- 
trast, invariance of f(v) with respect to O r is enforced 



by construction of the strain tensor y_ (a scalar in Sr), 
whereas that with respect to Q_ T is enforced by only al- 
lowing terms in f(v) that transform like a scalar under 
Q_ T , i.e., requiring that all the target space indices be 
contracted. 

As usual, the reference state, relative to which u is 
defined is taken to be in mechanical equilibrium, guar- 
anteeing that no terms linear in u appear. The first two 
terms of / are the standard elastic energy of an isotropic 
medium with A and \i the Lame coefficients |24| . We have 
included stabilizing nonlinear terms in the strain tensor 
u because we will eventually want to consider phase tran- 
sitions to an anisotropic state induced by a decrease in 
the shear modulus fi below a critical value. In the spirit 
of Landau theory of phase transitions [[Tcfl, at present, we 
view fi as a phenomenological parameter that is allowed 
to vary and even become negative. As we shall see in 
more detail in Sec. IV, the origin of a diminishing fi in 
liquid-crystal elastomers is the instability of the isotropic 
state toward the development of nematic liquid-crystal 
order characterized by the Maier-Saupe order parame- 
ter Q. In /, we have left out one third-order and three 
fourth-order terms permitted by symmetry, namely ones 
proportional to (Tra) 3 , (Tra) 4 , (Tru) 2 Tru 2 and TruTru 3 , 
respectively. Though these terms can easily be included, 
their effect is small for the nearly incompressible systems 
of most interest to us. 

Our primary interest is in the state with spontaneously 
broken rotational symmetry that is produced when /1 falls 
below a critical value. To describe this state and the 
transition to it, it is useful to decompose Uij into its scalar 
(in Sr) and symmetric-traceless parts: 



\5ijUkk 



where 



2 ®ij V>kk • 



(2.14) 



(2.15) 



Using Eq. ( 2.14 ) in Eq. ( 2.13 ) and keeping only the low- 
est order terms in Tru, we obtain the model free energy 
density that we will use in discussions of the anisotropic 
state and the transition to it: 



with 



/ = iB[Tra - (E/B)Ttu 2 } 2 + /1 (2.16) 



/1 = iATra 2 - CTru 3 + L>(Tru 2 ) 2 , (2-17) 



where A — 2/i, B = A + |/i is the bulk modulus, E = 
E' — C, D = D' — E 2 /(2B) and for simplicity we have 
dropped qualitatively inconsequential cubic and quartic 
terms in Tru. |43]| 



C. Anisotropic systems 

Often the reference state is a crystal that is invariant 
only under operations of some subgroup of O3. In this 
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case, there are additional combinations of the strain ten- 
sor that are invariant under the reduced set of symmetry 
operations of Sr, and the elastic energy is in general de- 



scribed, to harmonic order in 



in terms of a 4th rank 



elastic-constant tensor Cijki, with / = ^CijkiUijUki- 
will be particularly interested in uniaxial systems with 
axis along no, for which the general form of the elastic- 
constant tensor is (but see Sec. Ill) 



We 



Cijki = CinoiUojnoknoi + C 2 (n in jSll- +n ok n M 5® 



(2.18) 



where <$■?• = <5y — notTioj. The elastic energy in three 
dimensions with the z-axis chosen along no is 

/uili = ^ClU 2 zz + C 2 U ZZ {U XX + Uyy) + \Cz{U XX + Uyy) 2 

+Ci{ul x + u 2 yy + 2u%) + C 5 {u 2 xz + u 2 yz ). (2.19) 

The strain uy is still invariant under arbitrary rotation 
in St, so / un i is invariant under these rotations, as it 
must be. The reduced symmetry of the reference state 
introduces an asymmetry between the reference and tar- 
get spaces, and it is no longer so useful to introduced the 
alternative strain tensor v unless we wish to discuss ex- 
plicitly coupling between strain and another target-space 
tensor-field order parameter, such as the Maier-Saupe or- 
der parameter for a nematic. 

The elastic energy / uni of Eq. (2.19) is harmonic in the 
nonlinear strain uij. Higher order terms in Uij are, of 
course, permitted and are in fact necessary to preserve 
full rotational invariance in Sr, which is present (but 
hidden), if the uniaxial asymmetry arises as a result of 
the spontaneous symmetry breaking of an isotropic state, 
as happens in nematic elastomers, introduced in Sec. |TJ. 

In semi-soft elastomers [^8|,0 , the rotational invariance 
of the soft-elastomer isotropic state (discussed next) is 
only weakly broken. Any model describing these systems 
must introduce anisotropy in such a way that both the 
isotropic and the anisotropic soft phases are reproduced 
when the anisotropy is set to zero. The simplest such 
model can be constructed by adding an anisotropic term 



Ja 



-hnoiUijn j, 



(2.20) 



which breaks Or symmetry, to the free energy of Eq. 
(2.16). Here no is a vector in Sr that specifies the direc- 
tion of preferred alignment, and h is a field measuring the 
anisotropy strength. The properties of this model will be 
explored in a separate publication |34f|. 



III. STRAIN-ONLY MODEL OF NEMATIC 
ELASTOMERS 



of Eq.( 2.13| ), there can be a transition from an isotropic 
state with A ~ 5 to a uniaxial one with two rather than 
one distinct eigenvalues for A. This nematic-gel state is 
obtained from the isotropic one by stretching or com- 
pressing along some arbitrary direction in Sr specified 
by a unit vector no, which without loss of generality we 
take to be along the z axis. It is characterized by an 
anisotropic equilibrium right strain tensor u Q with prin- 
cipal axis along no. The transition to the nematic gel can 
thus be described completely in terms of the free energy 
f(u). Alternatively, the nematic gel can be characterized 
by an anisotropic equilibrium left strain tensor with 
anisotropy axis along some unit vector ni is St, and the 
transition to it can be described by f(v). The nematic 
gel breaks both Or and Ot symmetry. The description 
in terms of u displays explicitly the broken Or symmetry 
and that in terms of v the broken Ot symmetry of the 
nematic gel. The underlying order parameter, however, 
is the deformation tensor A, which exhibits both broken 
Or and Ot symmetry in the nematic gel. Even though 
the nematic gel breaks two symmetries, they are both 
broken at the same time, and there is only one transition 
from the isotropic phase to the nematic gel. As discussed 
in Sec. |II B| , f(u) and f(v) are identical functions of their 
arguments, and u and v develop nonzero anisotropic val- 
ues simultaneously. 

Though free energies expressed in terms of the strain u 
and v provide complete descriptions of the phase transi- 
tion to the nematic gel, it is important to remember that 
the full position function R(x) or equivalently the dis- 
placement u(x) is needed to describe all configurations 
of the gel. The tensors A, u, and v only provide infor- 
mation about long-wavelength distortions. A full statis- 
tical mechanical treatment of nematic gels requires the 
inclusion of curvature energies depending on the second 
derivative of R(x) into the elastic energy that appears in 
the partition function trace. This will be discussed more 
detail in a separate publication p6f . 

In this section, we will explore the properties of the 
spontaneously formed nematic gel described in terms of 
u and f(u). The description in terms of v is essentially 
equivalent. We will explicitly derive the soft elasticity of 
nema tic gels whereby the strain elastic constant C5 [Eq. 
(J2.19[) 1 vanishes identically (T^j and there is zero stress 
p5[ associated with appropriate strains up to a critical 
value perpendicular and parallel to no as long as other 
strains [ ^4|j45| ] are allowed to relax to their lowest en- 
ergy configurations. Our treatment provides a complete 
description of nematic gels and transitions to them with- 
out any reference to undelying nematic order. In the next 
section, we will consider nematic order and its coupling 
to strain and show that instabilities toward the devel- 
opment of nematic order drive the decrease in the shear 
modulus discussed in the preceding section. 



Under appropriate conditions, for example for suffi- 
ciently small shear modulus fi in the model free energy 
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A. Description in terms of Uij 

It is quite clear from the c ubic form of the elastic free 
energy, Eqs. ( 2.13 ) and ( [2.171 ), tna t when \x becomes suf- 
ficiently small, for finite C, there is a first-order tran- 
sition from an isotropic to a uniaxially distorted elastic 
state, which is very similar to the familiar isotropic-to- 
nematic transition |j], 3. We will consider this transition 
in more detail in Sec!" ]!!! G . In this subsection, we will 
investigate the resulting anisotropic elastic state, whose 
properties depend only on the existence of spontaneously 
formed anisotropy and not on any particular model of 
the isotropic-to-nematic transition. 

In the positive (negative) uniaxial state that results 
from such transition, the elastic material is stretched 
(compressed) along n in Sr and compressed (stretched) 
along directions perpendicular to no. This anisotropy 
axis can point in any direction ni in St- For the mo- 
ment, we will take ni = no, i.e., we do not rotate the 
sample after it has been stretched. In this case, the coor- 
dinates of its mass points in St are, therefore, described 
by 



of the strain tensor u from its new equilibrium value 
w Q , both measured in the coordinates x of the original 
isotropic state, Sr. It is, however, more common and 
convenient to describe these distortions in terms of dis- 
placements R'(x') = R(x) and strains u'(x') expressed 
as functions of the coordinates x' = Ro(x) of the new 
equilibrium stretched state: 



R'(x') 



u'(x') 



l' = 5(A' T A'-|) 



+ u + £u(x) 



(3.7a) 
(3.7b) 



where Ay = dB^/dx'j and rj^ = du'Jdx'j. Since 

^•3 = «=4% (3.8) 



dx' k dxj 



the strain deviation Su is directly proportional to u' and 
therefore proportional to the symmetrized strain (rf + 

r)' T )/2 when linearized: 



^ = 4oi'4 - 



(3.9) 



R (x) = x' = x + u = A X , 



(3.1) 



where the deformation tensor is spatially uniform and 
given by 



A 0ij = A 0± 6 ij + A \\n 0i noj, 

= A j_% + (A || - A ±_)n 0i n j, 



(3.2a) 
(3.2b) 



where 644 = c 
librium strain tensor is 



noinoj. The corresponding right equi- 



u 
=0 



|(AlA 



=0 



5)- 



(3.3) 



The anisotropy of the uniaxial state can be characterized 
by the anisotropy ratio 



A 2 
A Q|| 

A 2 ' 



Since 



AoA^ A 2 ± [|+(r-l)n n ], 



(3.4) 



(3.5) 



it is clear that the system is isotropic if r — 1 and only 
anisotropic if r 7^ 1. Both positive (r > 1) and nega- 
tive (r < 1) uniaxial anisotropies are possible, but we 
will focus mostly on positive uniaxial systems. In incom- 

2 - 1 A - A" 1 / 2 

0_L - A 0_L - A || , 



pressible systems, det A o = A ||A 



and r = Aqm. Many of the properties of the uniaxial 
phase depend critically on r. 

The goal of this subsection is to explore the elasticity 
of this spontaneously uniaxially-ordered gel. Distortions 
such a system can be described by deviations 



% = !-io = l(A T A-A A ) 



(3.6) 



One would normally expect the elastic free energy for 
strains u' about the new uniaxial state to have the form 



of Eq. (2.19), characterized by five independent elastic 



constants. However, as discussed in the Introduction, the 
fact that the uniaxial state arose via spontaneous sym- 
metry breaking of an isotropic state guarantees that the 
shear modulus C5 must vanish. We now demonstrate this 
explicitly. 

Since the original free energy is invariant under ro- 
tations O in Sr. the anisotropy direction no in Sr is 

IX 

arbitrary, and states characterized by strain Q^u^ O^ 1 
and «q must have the same bulk energy. This means that 
there is no bulk energy cost associated with a strain 



i'W = (A^)- 1 (2^2- 1 -i )A- 1 



(3.10) 



relative to the uniaxial state characterized by u since it 
describes a rotation in Sr, and is therefore a Goldstone 
mode of broken Or symmetry. For rotations through 9 
about the y-axis, 



O, 



cos sin 9 

1 
— sin 9 cos 9 



(3.11) 



Using this R inside Eq.( 3.1C ), we find that u'^iff) is a 
nontrivial strain even though it describes a pure rotation 
in Sr. Under this rotation, u'^0) only has components 
in the xz-plane (the plane of rotation 9) that are 



^=^(r-l) 



1 - cos 29 
r" 1 / 2 sin 29 



r- 1 / 2 sin 2(9 
r- 1 (l-cos26») 



(r - 1) ( 
2v^ 



(3.12a) 
(3.12b) 



8 



where r is the anisotropy ratio introduced in Eq. (3.4) 
and the final form is valid for small 9. Since, as just 
argued, the elastic free energy must be invariant under 
rotations in Sr, it cannot depend on the rotation angle 
9, and, therefore, there must be no energy cost associ- 
ated with an infinitesimal strain u' xz — u' zx . Similarly, 
invariance with respect to rotations about the x axis im- 
plies no energy cost associated with the strain u'. Thus, 
the shear elastic modulus C5 must identically vanish in 
a spontaneously uniaxial state, whose harmonic elastic 
energy, 



TV 
uni 



|Ci?4 2 2 + C 2 u' zz (u' xx + u: 
Ci{u\ 



) + hC 3 (u' xx +u') 



12 

XX 



yy 



2u' 



XI) I 



yy ' 
(3.13) 



is characterized by only four elastic constants. The su- 
perscript N in f^ ni is introduced to distin guish it from 
the standard uniaxial energy / uni [Eq.( 2.19| )] with C5 ^ 0. 
Because f^ ni contains only quadratic terms in strain ir- 
relative to the broken-symmetry uniaxial phase, as in 
similar systems J47| it is only invariant with respect to 
infinitesimal rotations in Sr and terms nonlinear in up - 
most be incorporated in order to encode the full Or in- 
variance flisf . 

There are striking experimental consequences of th e 
existen ce of zero-energy strains v! (9) given by Eqs . ( 3.10 ) , 
(3.12a) for arbitrary 9. Namely, if one of the components 
of strain u' xzl u' xx , or u' zz is imposed with the right sign, 
the other two components can (boundary conditions of 
the experiment permitting) adopt values to produce the 
zero-energy rotational strain of Eq. (3.12a). When r > 1, 
this relaxation is possible only for positive u' xx (extension 
perpendicular to the uniaxial direction no), negative u' zz 
(compression along no), and either positive or negative 
u' xz . For negative anisotropy systems (r < 1), the zero- 
energy strain relaxation is possible only for negative u' xx 
and positive u' zz . 

To illustrate this, consider first a sample with n p 
aligned along the z-axis, with r > 1. From Eq. ( 3.12a ), 
it follows that u' xx — (1/2) (r — l)sin 2 # is positive for 
a rotational strain when r > 1. Thus, we can only have 
soft elasticity for extensional strains along x, and we take 
u' xx > 0. If no relaxation of strain is allowed, this stretch 
would cost an energy proportional to u xx . If, however, 
strain relaxation is allowed, strains 



1 



VKJr - 1 - 2u 'xx) 



(3.14) 



The angle 9 specifies the direction of the induced uni- 
axial equilibrium stretch axis relative to the z-axis of a 
fixed coordinate system in Sr. In the current problem, 
this anisotropy axis is initially along z, and it rotates to- 
ward the re-axis as u' xx is increased until, at the critical 
strain u' xx = (r — l)/2, 9 = tt/2 and the anisotropy axis 
has been rotated to be along the ic-axis, as illustrated in 
Fig. H For strains u' xx larger that (r — l)/2, the sample 
will merely stretch along its new anisotropy axis along x 
with the additional strain Su' xx = u' xx — (r — l)/2. We 
can calculate the energy associated with this addi tional 
strain from the harmonic free energy of Eq. ( |3.13| ) pro- 
vided we remember that Su' xx is measured relative to the 
original reference system with the anisotropy axis along 
z rath er than x, i.e., if we remember to replace u' zz in Eq. 
O by u' zz = (A 2 JA*)5u' xx = Su' xx /r 



u' xx and u' yy 



T3|) by u' zz = (AQ ± /A" oll )du xx 
should be rescaled as well, but since we minimize over 
these quantities at fixed Tt zz , we do not have to explicitly 
consider these rescalings. Performing this minimization, 
we find 



Sf = 



0, 



2T 1 \ 



Ci 



Cr, 



2C4+2C3 



(*4 



if 8u' x: 
if 6uL 



<0; 
> 0. 

(3.16) 



Consequently, the stress df/dA' xx — A' xx df /du' xx for an 
ideal nematic gel is zero for u' xx < (r — l)/2 and grows 
linearly in 8u' xx for u' xx > (r — l)/2 as shown in Fig. 




X 



" XX XX 

FIG. 3. Schematic representation of the transformations of 
a rectangular piece of soft nematic gel subjected to a strain 
perpendicular to its anisotropy axis. If u' zz and u' xz are con- 
strained to be zero, there is no rotation of the anisotropy axis 
indicated by a double arrow as shown in the middle figure. 
If these quantities are allowed to relax, the anisotropy axis 
rotates to produce a state with the same energy as the ini- 
tial state as shown in the final figure. In the process, the 
original rectangle is transformed into a parallelogram sheared 
through and angle 8e, and its anisotropy axis represented by 
the double arrow is rotated through an angle 9t- 



convert the 



strain to a zero-energy rotation strain 



tensor (Goldstone mode) with rotation angle 



2u' 



(3.15) 



Thus, in an ideal system, there is no bulk energy cost 
associated with strains < u' xx < (r — l)/2. 
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A' 



FIG. 4. Stress a xx versus strain A xx for an ideal soft ne- 
matic elastomer stretched along a direction perpendicular to 
the direction of initial alignment. The stress is zero up to 
a critical strain = Beyond that, the stress initially 
grows linearly from zero. 

The angle 9t of rotation of the anisotropy axis in the 
target space St is not the same as the rotation angle 9 
in Sr. Indeed, since the energy is invariant with respect 
to rotations in St, 9t would be arbitrary, were it not for 
boundary conditions. A specific experimental geometry 
might for example demand that there be no change in 
the z coordinates of mass points as a function of x' . This 
is the situation depicted in Fig. |^. In this case, A' zx = 0, 
but A' xz ^ 0, and the xz submatrix of A^- takes the form 



A' 



A' A' 

XX XZ 

A' 



(3.17) 



This deformation tensor corresponds to a zero-energy ro- 
tation in Sr described by a rotation matrix Q_ R [Eq. 



(3.11)] provided there exists a rotation matrix 



cos ot sm &t 



(3.18) 



sm t>T cos Ot 

in the xz-plane of St such that the strain A = A'A 
relative to the original reference system is OjAqO" 1 , or 



A' = £ T A |2~ 1 A ~ 1 , 



(3.19) 



an expression that is fully consist ent w ith the form of 
strain Goldstone mode Eq. 3.12a. A straightfor- 

ward calculation yields 



A' 



cos sm Vt 



<r sm cos Ot, 



(3.20) 



which, upon imp osition of the boundary condition of 
A' zx = 0, Eq.(p%, gives 



sin 9t = 



r sm 



1 + (r - 1) sin 2 
/ V 1 A 7 ^ 



(3.21a) 
(3.21b) 



where we used sin 9 = 2u' xx /(r — 1), Eq. ( 3.15 ) and 
(A^) 2 — 1 = 2u' xx to obtain the final form. This is ex- 
actly the same result obtained via a direct minimization 
of the neoclassical rubber energy in the incompressible 
limit [0. 



The angle 9t is simply the angle that the uniaxial 
anisotropy axis makes with the z axis in the target space. 
This axis is the principal axis of AA T (or equivalently of 
v) . A direct calculation of A'A Ag A' T yields 



AA 



^2 / cos 2 6>t + r sin 9t (r — 1) sin 9t cos 9t 
0J - ' (r — 1) sin 9 T cos 9 T r cos 2 9 T + sin 2 9t, 



(3.22) 



which is nothing more than AoAo rotated though 9t 
so that the principal axis is along ni = (sin 9t, cos O t)- 



Under the transformation defined by A' of Eq. (3.17), a 
rectangle will be transformed into a parallelogram with 
two edges parallel to the x-axis and two edges making an 
angle 



tan 



A' 



tan 



(r - 1) tan 9 T 
r + tan 2 9t 



(3.23) 



with the vertical. Both 9e and 9t are indicated in Fig. 
||. Note that the argument of the inverse tangent in the 
expression for 9e rises from zero, reaches a maximum, 
and then returns to zero as 9t passes from zero to 7r/2, 
the maximum angle of rotation of the anisotropy axis. 
In the process, A passes from having extension Aqii along 
the z-axis and Aoj_ along the x-axis to having extension 
Aoj_ along the z-axis and Aqii along x. Thus, the orig- 
inal rectangle is distorted to a parallelogram that first 
becomes more slanted as strain increases, reaches a max- 
imum slant, and then becomes less slanted until it finally 
reaches a rectangular form that is precisely the original 
rectangle rotated by 7r/2 before stretching further along 

the x axis. 

It is notable that the expression, Eq. (|3.21b|) , for sm Ot 
is independent of the detailed form of the elastic energy 
and that it is characterized only by the level of anisotropy 
r of the initial uniaxial state and by the deformation A' xx 
applied to it. The advantage of our approach is that 
it makes it clear that the phenomenon of soft response, 
summarized by Fig. ||, i.e., zero stress for a range of lon- 
gitudinal strain, applied perpendicular to the uniaxial di- 
rection, follows entirely from general symmetry principles 
of breaking of rotational invariance of the reference state. 
It is a property of any nematic gel formed spontaneously 
from an isotropic gel, and is therefore independent of the 
details of the microscopic model of the gel and the mech- 
anism that drives the uniaxial in stabilit y, as confirmed 
by the generic form for 9t in Eq. ( 3.21b ). 

There is spectacular experimental evidence |Q for 
-(A^), Eq. ( 3.21b| ). Likewise, experiments confirm the 



GL prediction of softness (vanishing stress up to a criti- 
cal value of strain, Fig. ^) , in the limit in which the gel is 
crosslinked in the isotropic phase. In contrast, chemically 
identical networks crosslinked in the nematic phase ex- 
hibit a plateau in the stress-strain curve that approaches 
zero as the degree of nematic order during crosslinking is 
decreased. See Ref. [^0| for a presentation and discussion 
of these results. 
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It is clear that for r > 1 an imposed shear strain u' xz 
or a compressional strain u' zz < can be converted to 
zero-energy rotational strain just as for the case u' xx > 
just considered. The energy of a shear strain is zero for 
\ u 'xz\ < { r ~~ l)/(4-\A")> and that of a compressional strain 
is zero for \u' zz \ < (r — l)/(2r). While softness with re- 
spect to small shears below a critical frequency has been 
observed fl48|| , to our knowledge, our above prediction of 
softness over a finite range of strain for the geometries 
with imposed u' xz ^ or u' zz < has not been tested. 

The above discussion a nd th e structure of the zero- 
energy strain Eq. fl3.12a ), imply that our argu- 
ments for soft elasticity go through equally well for the 
negative uniaxial anisotropy elastomers, r < 1, but with 
reversed signs of the imposed zero-energy strains. Thus, 
in negative uniaxial elastomers, there is no stress asso- 
ciated with unconstrained compression along x, u' xx < 
and extension along z, u' zz > 0. 



B. Description in terms of Vij 

In the discussion just presented, the nematic state and 
its soft elasticity were described in ter ms of the right 



strain tensor mj. As demonstrated in Sec.|[IB|, this state 



could equally well be described in terms of the left strain 
tensor Uy, which has an equilibrium value 



Eo= 2 



1(4^-1) 



(3.24) 



that is identical to u Q in the basis defined by no. The de- 
viations from the spontaneously anisotropic equilibrium 

state A can be described by 
=o J 

^ = I-io = 5(M T -A A^) (3.25) 

as well as Su. The deviations Su and Sv, however, have 
different relations to the displacement gradient tensor rf 
relat ive t o the anistropic equilibrium state. As we saw in 
Eq. ( |3.9| ), Su is linearly proportional to u' and at linear 
order depends only on the symmetrized part of rf . Sv 
on the other hand is not proportional to v' defined by 
Eq. (2.7) with A replaced by A', and to linear order, it 



depends on both t he sym metr ic an d anti-symmetric parts 
of r/' . Using Eqs. ([3.25]) and ([2.2]), we easily derive 



+l( r ~ l)(«oi«ofe7?jfe + Vik n oknoj 
+v'ik n okn ir)' jl ), (3.26) 

in terms of rf , with r the previously defined anisotropy 
ratio [Eq. (|3~4|]. Using the decomposition 

(3.27) 



rf = rf + rf 

= =S =A 



IT I I 

2 = R s ~ R A ' 



where if Sij = rf Sji and rf M 



i , ^^Aji are ' respectively, the 
symmetric and anti-symmetric parts of rf^ and defining 
a rotation angle 



i du 'k 



l^ijkVAkj 



(3.28) 



we obtain 



A 0± 5 Vij 



+ ±(r - l)[n 0i n kV'sjk + v'sik n okn j] (3.29) 
-§(r - l)[n 0i (n x Sl)j + (n x Sl)in j} 

to linear order in rf^. We thus explicitly demonstrate 
that Svij is a function of both the symmetric and anti- 
symmetric parts of rjjj. At first pass, this observation ap- 
pears to co ntrad ict t he fa cts that the linearized form of 
Suij, Eqs. ( 3.7b ) and (3.S), does not depend on rf Ai p and 
that the harmonic free energy of the nematic phase must 
have exactly the same form whether expressed in terms of 
Su^ or Svij . The dilemma is resolved by noting that only 
the soft components of the strain, 5v xz and Sv yz depend 
on rf Ai p and that, to harmonic order, f^ ni (v) is guaran- 
teed by Or invariance (which leads to vanishing of C5) to 
be independent of such strains. Consequently, consistent 
with the expectations, in both tty and descriptions 
no anti-symmetric part of the deformation tensor rf Ai A 
appears. 



C. Isotropic-to- uniaxial transition 

In a transition from the isotropic to the uniaxial state, 
the strain develops a nonvanishing anisotropic compo- 
nent. We can describe this transition equivalently in 
terms of u and f(u) or in terms of v and f(v). To 
be concrete, we will use u-description here. Since the 
isotropic part of the strain is insensitive to anisotropy, 
the appropriate order parameter for the transition is 
t he sy mmetric-traceless component of the strain, u [Eq. 
( |2. 15 )], which is identical in form to the symmetric- 
traceless order-parameter tensor Qij of a nematic liquid- 
crystalline phase Thus, in mean- field theory, the 
transition from the isotropic to the uniaxial state is iden- 
tical to the isotropic-to-nematic transition, whose prop- 
erties have been exhaustively studied [QJ|] . (In Appendix 
[A|, we review the formal properties of this transition that 
are relevant to the current discussion.) To see this in 
more detail, we ca n integrate out the "massive" Tru from 
f(u) in Eq. ( 2.16 ), to obtain an effective theory in terms 
of u alone. This operation yields 



E 

Tru = — Tru , 
= B — 



(3.30) 



and t he effective free energy reduces to fi(u) of Eq. 
( |2.17 ). Because of the presence of a cubic invariant, 
the free energy fi(u) exhibits a first-order transition at 
A = A c = C 2 /{12D) to a state with 



u% = V(r°~° 



I A. 

3 °y > ' 



where ip satisfies the equation of state 



(3.31) 
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A - Gib + %Di\) 2 = 0. 

The total strain in the distorted state is thus 

2E 



4 = 4 



9B^ 



(3.32) 



(3.33) 



This corresponds to a stretched state with target space 
positions Ro = A X characterized by a deformation ten- 



A = Jl + 2u? = 



A ± 
A ± 
A || 



(3.34) 



with 



A ± = ^1 
A o|| = f 1 



AE 
9B 



ib 



ib 2 -\ — ib 

95 3^ 



1/2 



1/2 



This form for A Q preserves the volume up to order ip 2 . 
The order parameter ib is a direct measure of the spon- 
taneous stretch anisotropy of the nematic state, with 



ib 



2"( A 0|| 



A ox) 



h A o±( r 



1) 



(3.36) 



where r is defined in Eq. (3.4. 

The elastic free energy fi(u) can now be expanded in 
powers of Su = u — w Q and reexpressed in terms of vf , de- 
fined in Eq. (|3.9|). Ward identities, imposed by the rota- 
tional Or invariance, guarantee that terms proportional 
to (5u xz ) 2 and (Su yz ) 2 vanish. The h armo nic elastic en- 
ergy is, therefore, given by f^ ni [Eq. ( 3.13 )] with 



established uniaxial order, i.e., the uniaxial state be- 
comes unstable relative to a biaxial state with different 
equilibrium strains in all three directions. A biaxial ne- 
matic gel is softer than a unaxial one ]5M , and, as we 
will show here, it has no nonvanishing shear modulus in 
three dimensions. The order parameter for the unaxial- 
to-biaxial transition is the two-dimensional symmetric- 
traceless tensor obtained by projecting u onto the xy- 
plane. Since there are no cubic invariants of a two- 
dimensional symmetric traceless tensor, the transition 
from the uniaxial to the biaxial state is generically a con- 
tinuous transition in the xy universality class. 

The biaxial phase is characterized by a Cauchy defor- 
mation tensor with three independent components: 



4> 



Aqi 
A 02 
A 03 



(3.38) 



(3.35) and the corresponding equilibrium strain tensor given by 



uoi 
■g = I u 02 
«„ 3 



(3.39) 



with u 0a = (Ag Q - l)/2, a = 1,2,3. The additional 
broken rotational symmetry of the biaxial relative to the 
uniaxial phase causes more shear elastic moduli to van- 
ish. As in uniaxial gels, strai ns of the form of Eq. ( 3.10| ) 
[with m q given by Eq. (|||)] that arise from arbitrary 
three-dimensional rotations Or in Sr cost no energy. For 
simplicity, we consider only the most general infinitesi- 
mal rotation matrix, which can be expressed in terms of 
rotation angles 9 X , 6 yi and 6 Zl respectively about the x, 
y, and z axes as 



Cx - 
C 2 = 
C 3 = 



AE 

11 ' ' ~ 3^ 



4 2 
-A+-Cib 
3 3 v 



A 



oil 



AE \ ( 2E \ 2 , 1 



B (l + —ib) -lA--Cib 
{ 2>B ) 3 3 v 



A 0_L 



C 4 = -ObA* x , 



A 2 \2 
A 0|| A 0_L 



(3.37) 



with C5 = 0, as anticipated in our early discussion of 
the generic, symme try-d ictated form of the clastic free 
energy. From Eq. ( 3.32 ) we note that ib has the same 
sign as C, ensuring that C4 ~ Cip is always positive. 



D. Biaxial nematic 



It is clear from the form of f^x(u'), Eq. ( 3.13| ) that 
if C4 is driven negative Q] the uniaxial state becomes 
unstable to strains in the xy-plane perpendicular to the 



1 -6 



Or 



1 -o x 

0* 1 



(3.40) 



The symmetric zero-mo de st rain tensor yf calculated 



from this Or and Eq. (|3.10| ) has components that to 
linear order in the infinitesimal angles of rotation are 
u' — u' — u' — and 



= AqiAq2(u 01 - u 02 )9 z 
= -Aq/A^Cwqi - u 03 )9 y 
= Aq2 1 A^(u 02 - u 03 )6 x . 



(3.41a) 



Since the underlying Or invariance demands that there 
can be no energy cost associated with such zero-mode 
strains, the elastic energy cannot depend on the shear 
strains u' xy , u' xz , or u' to harmonic order. The hallmark 
property of solid that it can support a static shear stress 
is therefore lost in a spontaneously biaxial solid. The 
biaxial nematic is an anisotropic tethered fluid [|l] |53| . 
The harmonic elastic energy of a biaxial gel, therefore, 
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depends only on the compression/extensional strains u' x . 



u yy and u' zz and has the form 



/biax 2 ^ , ^ 



a(3 u aa u t3p- 



(3.42) 



There are in general six independent components of B a p. 
As in uniaxial gels, there is soft compressional and exten- 
sional elasticity in biaxial gels with vanishing stress up 
to critical values of the strain. We will not treat these 
properties in detail here. 



IV. NEMATIC GELS: STRAIN AND 
ORIENTATIONAL ORDER 

In this section we extend our formulation of the model 
of anisotropic gels to include both the elastic and ori- 
entational (nematic) degrees of freedom. We first con- 
sider a "soft-spin" theory in which orientational order is 
described by symmetric-traceless nematic order param- 
eter Qij, which has both uniaxial and biaxial compo- 
nents. This theory, which can describe both isotropic 
and anisotropic phases of gels and the transitions be- 
tween them, takes explicit account of the coupling be- 
tween strain and Qij. It can be viewed as a theory 
in which the familiar isotropic-nematic transition char- 
acterized by ordering of Qij induces elastic distortion. 
Guided by the underlying rotational symmetry of the 
nematic gel, we then develop a complementary "hard- 
spin" model of nematic gels valid deep in the nemati- 
cally ordered phase. This theory is formulated in terms 
of the strain and the nematic director n alone, with all 
"massive" modes (e.g., magnitude of the uniaxial order S 
and biaxial fluctuations) integrated out. A common fea- 
ture of these complementary models is their invariancc 
with respect to global simultaneous rotations of strain 
and nematic order. This invariance leads to gauge-like 
couplings between strain and nematic order, whose har- 
monic limit reduce to those derived by Olmsted jl4| fol- 
lowing de Gennes ||. However, our expression in terms 
of Vij and n for the globally invariant energy deep in the 
nematic phase is new. 



A. Simple model of the IN transition 

In the preceding section, we investigated a model in 
which an isotropic elastic medium undergoes a sponta- 
neous anisotropic distortion triggered by the fall of its 
shear modulus below a critical value. In liquid-crystal 
elastomers, the reduction in the shear modulus and the 
elastic distortion that it leads to are actually driven by 
the underlying isotropic-nematic transition of the meso- 
genic component of the gel, that is orientational order- 
ing of e.g., side-chain or main-chain nematogens. It is, 
therefore, of some interest to develop a model in which 
the orientational order parameter Qij explicitly appears. 



A generic model free energy density for such a model of 
a liquid-crystal gel will consist of an isotropic elastic term 
fei(u), a term fq{Q) for nematic orientational order, and 

a nemato-elastic term /cfe Q) that couples strain to the 
nematic order parameter Qy : 



/cl-Q — /ol + fo + fc- 



(4.1) 



For simplicity, we can take f c \ to be the elastic energy 
/ of Eq. Q2 . 1 3[ ) with only quadratic-order terms in Uij 
(or, equivalently, Vij), and near the IN transition, we can 
choose the usual Landau-de-Gennes form for 



fc 



Q ~ 2'Q 



TiQ - w 3 TrQ + ^(TrQ 2 ) 



2\2 



(4.2) 



Terms in gradients of Q should also be included, but 
they do not affect the present mean-field discussion, and 
we will therefore ignore them here. The most gen- 
eral local energy coupling strain to Qij can be con- 
structed from products of terms invariant under arbi- 
trary rotations in both Sr and St, whose general form 
is Ti[g ni Q mi ...g 1 pQ mp }. Note that these terms involve 

couplings between (rather than ) and Qij , because 
Qij exists in St , and like transforms like a tensor un- 
der rotations in St but like a scalar under Or rotations 
in Sr. To keep our discussion simple, we will for the 
moment consider a simple form for fc' 



fc = -sTruTiQ - 2tTrvQ, 



(4.3) 



where V; 



of 



^SijVkk is the symmetric-traceless part 



and where we used the fact that Tru = Trv. 



This energy captures the important qualitative features 
of strain-orientational coupling, namely that the devel- 
opment of orientational order will drive an anisotropic 
distortion and a smaller change in volume. 

An elastic energy f£ that is a function of strain alone 
can be obtained by integrating Q out of the total free 



energy of Eq. (4.1). The leading order correction of this 
operation to f c \ is — 2(< 2 /r^)Tru 2 . Thus /"[ has exactly 
the same form as Eq. (2.13), with \i replaced by // = 
/i — (2t 2 /rg). Clearly /? decreases and passes through 
zero as r' decreases and the IN transition is approached 
from the isotropic phase. Thus, the decrease in fj, in 
the models of Sec. Ill arises from instabilities toward the 



development of nematic orientational order. 

To treat the effects of strain-orientational coupling af- 
ter the transition to the nematic state occurs, it is useful 
to recast / in a slightly different form: 



/d-Q = HTr i -( s /5)Tr£ 2 ] 2 
+IiTr{g-(t/n)Q} 2 + f Q 



where 



2\2 



f Q = ±r Q TrQ 2 - w 3 TrQ 3 + Wi (TiQ 2 ) 



(4.4) 



(4.5) 
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with r Q = r' Q - 2(t 2 /fi) and w 4 = w' 4 - (s 2 /2B). This 
free energy leads to the equations of state 



df 

df_ 

dv i:j 

df 

dQ.,3 



= B[Ti^- (s/B)TrQ 2 } = 

= uliiij - (t/n)Qij} = (4.6) 

-2t[v l3 - {t/n)Qij] = 0. (4.7) 



The uniaxial solutions to these equations are given by: 
s 



Tru 



59. 
> j 



-TrQ^, 
B =o' 

Qij 7 



Q° tj = S(n 0l n 0j - ISij), 
with S satisfying 



(4.8a) 
(4.8b) 
(4.8c) 



r Q S - w 3 


S 2 + ^w 4 S 3 = 0. 


(4.9) 


Using these uniaxial 


solutions Eq. (fD 


|) in the new 


stretched state, we find 




Ag_L = 1 + 


9B 3/i 




Ao|| = l + 


As s 2 + 4t s 

9B b + 3v 


(4.10) 



Note that — Aqj_ = 2(t/fi)S is linear in the nematic 
order parameter S. 

As discussed in Appendix |a|, fluctuations away from 
the equilibrium state are conveniently treated with 
the introduction of a complete set of five orthonormal 
symmetric-traceless matrices satisfying I^lfi = S af3 
that allow us to expand Qij and Vij as Qij — X)q=o Q^Ifj 

and = ^2a=o Va ^ij- Expressions for Q a and v a in 
terms of Qij and uy, respectively are given in Eqs. ( |A2| ). 
In particular, Qq = ^2/3S. In terms of these variables, 
we have to harmonic order in Svij and 5Qij 

Sf*-Q = ±B[TrSg-{4,s/3B)S5S} 2 
+H[5v - (t/fi)^/2j3SS} 2 

4 

+/i^[« Q - (t/fl)Q a ] 2 
a=l 

+ ±A 1 {6S) 2 + ±A 2 [Q 2 + Q 2 }, (4.11) 

where A\ and A 2 are given in Eqs. (|A6| ). Rotational 
invariance of /q guarantees that terms Q 2 ~ Q xz and 
Q\ ~ do not appear in the nematic state. We can 
integrate out the "massive" longitudinal mode 8S and 
biaxial modes Q\ and Q2 to obtain 



Sf v = \Bi(8v zz ) 2 + B 2 Sv zz (Sv xx + Svyy) 

+ k~B 3 (5v xx + Sv yy ) 2 (4.12) 

+B 4 (Sv 2 xx + Sv 2 yy + 25vl y ), 

+2fi ([6v xz - (t/fi)Q xz } 2 + [Sv yz - (t/fi)Q yz } 2 ) , 

where the coefficients B a are evaluated in Appendix 
This free energy is manifestly invariant under arbitrary 
rotations in Sr because it is a function of the strain 
only. However, its invariance in the St is restricted 
to infinitesimal rotations in Ot because we only used 
the harmonic free energy to integrate over "massive" 
modes. Because underlying Ot invariance of the ne- 
matic state forbids "massive" terms in Q xz and Q yz , in- 
tegration over them also eliminates strains v xz and v yz 
from the resulting elastic free energy, which, as antici- 
pated takes the form identical to that in Eq. ( |3 . 1 3| ) . Such 
symmetry-enforced vanishing of an elastic constant (here 
C5) is mathematically closely related to the well-known 
Anderson-Higgs mechanism in gauge theories. 

The terms involving Q xz and Q yz are interesting be- 
cause they determine the energy cost of rotating the di- 
rector away from the direction of uniaxial stretch. When 
we convert to the strain variables of the stretched state 



using Eq. ( 3.29 ) and the expressions, Eq. ( 4.10 ), for A 



VO-L 



and A || in terms of S, we obtain to lowest order 
(5n = n — no 

fT = l C lV'izz + ClVSzziVSxx + VSyy) 

+ h C 3(VSxx + v'syyf + C 4VSxx + VSyy + 2 VSxy) 



E Wsaz-f3(Sn a -v'Aza)f, 



-J J 



where 







(r + 1) 



(4.13) 



(4.14) 



and elastic constants C a are related to the constants 
B a via C\ = A i|Si, C 2 = A 2 l ^A 2 t± B 2 , C3 = Aq ± B 3 , 

C" 4 = A* ± S 4 , and fi' = (r + 1) 2 A% ± ^. The form of this 
energy is in fact the most general one, and we will derive 
it again in the next subsection after we have derived its 
nonlinear generalization. It is exactly the form obtained 
by Olmsted H] following de Gennes || and Bladon, Ter- 
entjev, and Warner |2^] . It shows clearly how the director 
can relax locally to Sn a — i]' Aza + P~ 1 i]s a ,z to eliminate 
any dependence of the free energy on rj' Saz , i.e., to make 
C 5 = 0. 



B. Theory with strain and director 

We have just seen how the development of nematic 
order characterized by Qij leads to a stretched nematic 
elastomer with a soft elasticity. The formulation in terms 
of Qij is well suited to a description of the transition 
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from the isotropic to the nematic state. Deep in the ne- 
matic phase, the theory that best captures the effects of 
long- wavelength strains and variations in the direction of 
nematic order is one expressed in terms of strain and the 
nematic director n only, i.e., one in which fluctuations in 
S and in the biaxial part of Qij are integrated out. This 
theory, like others we have discussed must be invariant 
under both rotations in Sr and under simultaneous ro- 
tations of Vij and n in 5r- 

To construct a fully rotationally invariant theory deep 
in the nematic phase, it is convenient to introduce a lo- 
cal coordinate system defined by the orthonormal triad 
{e^e^e 3 = n} consisting of the local director n and 
two vectors e 1 and e 2 perpendicular to n. These vectors 
satisfy 



E 

0=1,2 



e°e a 



(4.15a) 
(4.15b) 



In what follows, we will adopt a notation in which Greek 
indices /1 and v will run over 1 to 3, and Roman in- 
dices a and b will run from 1 to 2, i.e., over the subspace 
transverse to n. The left strain tensor v can always be 
expressed in terms of its components in this basis: 



where 



and 



V \\± ~ n i v ij e j 



,.ab 



% 3 e.y 



(4.16a) 
(4.16b) 

(4.17) 

(4.18a) 
(4.18b) 

(4.18c) 



The components v >iV are invariant under rotations in Sr 
because Vij is invariant under Or by construction. They 
are also invariant under simultaneous rotations of both R, 
and the triad {e M } in St, i-e., they maintain their same 
numerical value, under simultaneous rotations of and 
the basis {e^e^n}. 

A gel whose anisotropic state forms via spontaneous 
symmetry breaking from the isotropic phase has no pre- 
ferred or imposed directions, and the elastic free energy 
will depend only on v^, , and v^ ± . Furthermore, this 
free energy cannot depend on the arbitrary choice of the 
vectors e 1 and e 2 in the plane perpendicular to n, and 
it will be a function only of v^ v in the combinations v\ \ , 
u5°, Wm^Um^, and v'^v'f 1 . Since linear terms proportional 
to D|| and are present in the anisotropic phase, it will 
be characterized by a nonvanishing equilibrium strain Vq 
with components Vq\\ and Vq^_. If the equilibrium director 
is no , then such a uniaxially distorted state is character- 
ized by the equilibrium strain 



where 



( A o|| 



VQ±(Sij 
2 L v 1 f 'J K^-o) 1 

l)/2, v ± 



Voij = v \\n in j 



noiUoj), 



(4.19) 



(A 2 



0-L 



l)/2 and 



Gy(n) = A^niUj + Ao ± (% - riirij). 



(4.20) 



Away from equilibrium, the free energy can be ex- 
panded in the deviations 



(4.21) 



of the strain from its equilibrium value. To harmonic 
order in these deviations, we have 



Sfw = \c lW l + C 2Wll wT + ±c 3 «°) 2 

+C 4 wfwf + C 5 w^ l± , 
where the rescaled invariant strains are 

-2/ 



(4.22) 



w \\ = A o\\( v \\ 



a b 



, -2 



a b 



- A 0± 
ft _I_ = A o ( v \\±. 



vo\\) 



V, 



) 

o||x) 



0-L 
V? 



(4.23) 



-r2 



with Aq = (Aqm + Aq_|_)/2. The prefactors in these defi- 
nitions have been chosen to produce, as we shall see, the 
simplest linearized version of the harmonic free energy. 
Even though this free energy is only an expansion about a 
uniaxially-distorted equilibrium state, unlike the free en- 
ergies [Eq.(3.13)] expressed in terms of strain alone that 
we considered in Sec. Ill, it is completely invariant with 
respect to arbitrary rotations in both Sr and St- 

The equilibrium values Vq\\ and Vq± do not depend on 
the particular direction of n. If {e M } = {e^}, then v^ v — 
e oi v oij e oj- Away from equilibrium defined by no, My = 
VQij + Svjj where voij is given by Eq. ( [1.19| ) and Svij by 
Eq. ( gj§ . Thus 

Sv^ = (efej 

and we have 



3 0i e 0j>0ij 



(4.24) 



w \\ = - |(1 - (l/r))[l - (n • n ) 2 ] 



w a j_ a 



#«y + i(r-l)[l-(n-no) s ] 
pb 6 Q± v' -=rf + rt 

*3 ' L J txx 1 fyy 

\\± w \\± = A i n i n k 5 G j j L 8v lj 5v k i 

+2/3A 2 (n • n )niSvij5f k riok 
+/3 2 (n-n ) 2 [l-(n-n ) 2 ] 



(4.25a) 
(4.25b) 



1 



Wszl - P( Sn l - VAlz)Y 



(4.25c) 



lt ab n ,,ab 



+ (r - l)noinojS ik S lj Svki 
+ i(r-l) 2 [l-(n.n ) 2 ] 
SikSjh'sijVski = € x + Vyy + 2r@x V , (4-25d) 
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where (5 is defined in Eq. (4.14). The energy 8f w is char- 
acterized by the five elastic constants C a and the stretch- 
ing ratio r, which has the same value in every one of the 
nonlinear strains. 

Alter nativ e but equivalent expressions for the strains 
in Eq. ( 4.23 ) are useful and elegant. The components of 
the equilibrium strains Vq U have the same value if the 
basis {e^} is transformed to the b asis {e''} provided the 
director no in GV,(no) [Eq. (4.20)] is transformed to n. 
Thus we have 



J 



(e&G^CnoK,-^) 



0j 

= \{etG ij {n)e)-5n, 



and from Eqs. (p.7|) , (|3.S| ), and ( |4.17|) 

^ = i[ e fA^G ife (no)A^er-^] 
^From this we obtain 



where 



X = 3 (A'g(no)A - g(n)) , 



(4.26) 
(4.27) 
(4.28) 
(4.29) 



and finally 



v \\ 


— A Q || riiVijUj 




w a j_ a 






w u w u 


4 . 

= A riiJifc^Vijyfcj 




wfwf 


= A ol S i~k$fiVijVki. 


(4.30) 



C. Crosslinking in the Nematic Phase 

If an elastomer is crosslinked in the nematic rather 
than the isotropic phase, the memory of the anisotropy 
of the state, with a uniaxial direction n , at the time of 
crosslinking is locked in, and full Or invariance of Sr is 
reduced down to D^h symmetry. If coupling to nematic 
order is turned off, the system will be cha racte rized by a 
uniaxial elastic energy of the form of Eq. ( 2.19|) with five 
elastic constants in general. (Turning off this coupling 
is not as unphysical as it may seem. This is precisely 
what is done in treatments of plastic crystals consisting of 
anisotropic molecules such as N%-) This part of the elastic 
energy is a function of Uij and is invariant under rotations 
in St- It is also invariant under simultaneous rotations 
of no and x in Sr and under operations on x in Dooh at 
fixed no. Couplings to the nematic order parameter Q 
must be invariant under simultaneous rotations of R and 
Q in St and under simultaneous rotations of x and no in 
Sr. The simplest couplings linear in Q are of the form 



where hij — hn Qi noj and, as before, AA T = 5 + 2v. The 
first term in this energy reduces to —hriQiQijnoj and fa- 
vors alignment of principle axes of Q along no in the 
absence of deformation, when the deformation tensor A 
is the unit tensor. 



The generalization of Eq. (4.22) to systems crosslinked 
in the nematic phase is fairly complicated. It cannot 
be expressed in terms of the strain Svij alone; it can 
only be expressed in terms of the more fundamental non- 
symmetrized strains r\ij. However, the major effect of 
crosslinking in the nematic phase is to to make no a pre- 
ferred direction with an energy cost to rotate away from 
that direction, which can b e des cribed by the addition of 
a term — h(no ■ n) 2 to Eq. ( 4.22 ) to lowest order in rjij. 



V. NEOCLASSICAL THEORY OF ELASTOMERS 

So far we have described liquid-crystal gels in terms of 
nonlinear strains, rotationally invariant in either Sr or 
in St, relative to some equilibrium reference state, and 
we have focussed on those properties that result from 
the spontaneous broken rotational symmetry of the ne- 
matic state. We have treated the elastic constants in 
our model free energy as phenomenological parameters 
to be determined experimentally. To date, experimen- 
tal realization of liquid-crystalline elastomers are cross- 
linked liquid-crystalline polymers. They are rubbers with 
orientational degrees of freedom of a liquid crystal, and 
their elastic properties over a very wide range of strains 
can be described quantitatively by a generalization of the 
classic theory of rubber elasticity j2lj. This is a semi- 
microscopic theory in which the origin of shear moduli is 
the reduction of entropy arising from constraining poly- 
mers to pass through cross-linking points. In this section, 
we will show that this theory, when expressed in terms 
of non-linear strains, is equivalent to those discussed in 
preceding sections of this article. 

In the simplest version of the neoclassical theory, poly- 
mer segments between crosslinks are viewed as indepen- 
dent random-coil polymers of length L. In the anisotropic 
environment induced by the nematic order, the effective 
step lengths parallel and perpendicular to the direction 
of nematic order are different, and mean-square end-to- 
end displacement is characterized by an anisotropic step- 
length tensor, 



I = ig, 



(5.1) 



where I is a length and g is a unitless tensor, reflecting 

system anisotropy, whose form will be discussed in dif- 
ferent contexts below. The probability that the two ends 
of a single chain are separated by R is 



P(R) = 



fg = -TrA h A r Q - 2/3TrvQ, 



(4.31) 



detr 



(2ttL/3) 3 



exp 



(5.2) 
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The free energy per chain is /chain = — TTnP(R). Now 
assume that the separation R was produced by an affinc 
transformation from some initial state with separation 



Ro such that iL 



A r ijRoj, where A is the deforma- 



tion tensor relative to the initial state. (Later we will 
introduce a new reference state and use the symbol A to 
denote deformations relative to that state.) The free en- 
ergy per chain of the entire elastomer is then /chain (Ro) 
averaged over all separations Ro of the initial state, which 
we assume consists of random- walk chain segments char- 
acterized by a step-length tensor i n = log and a proba- 



bility distribution given by Eq. (5.2) with I replaced by 



lq. The initial state may be viewed as the state at the 
time of crosslinking. Thus, if the system is crosslinkcd 
in the isotropic state, Z Q will be an isotropic tensor; if it 
is crosslinked in the nematic state at some temperature 
T, the degree of anisotropy of l Q will reflect the degree 
of nematic order at that temperature. The free energy 
density relative to the initial state is thus 



/, 



ch 



±nT(TrA I A 1 1 

2 *. =r=0= 



T,-l 



lndety -1 ), 



(5.3) 



where n is the volume density of chain segments. This 
purely entropic free energy, whose ground state is the 
collapsed state with A = 0, cannot alone provide a 
complete description of the elastic properties of an elas- 
tomer. It must be supplemented with some treatment 
of the short-range enthalpic forces that prevent collapse 
to infinite density. Merely imposing the incomprcssibility 
constraint, det A = 1, is sufficient to provide a very good 
description of dense nearly, incompressible systems. We 
will take a phenomenological approach in which there is 
an energy cost, measured by a compression modulus B r , 
arising from deviations of det A from 1 : 



/b = ! -B r (det A - If 



(5.4) 



Our complete neoclassical energy density is thus / = 
/ch + /b • 

An important feature of this model is that it depends 
on A^ only via the combination A^g A? because the de- 
terminant of a product of tensors is the product of the de- 
terminants. Thus, it is convenient to analyze this model 
in terms of A = A g 1 ^ 2 , the strain tensor relative to the 



isotropic state obtained by rescaling lengths via g ' . |5J] 
Our model is thus 



/ = |nT(Z //)TrM T |- 1 + \nT Xng^g 1 
+|S r [(detM T /det£ o ) 1 / 2 n2 



1] S 



(5.5) 



where n is the volume density of chain segments. 

We will now analyze two version of this model: one 
appropriate to the description of the IN transition and 
one appropriate to systems deep in the nematic phase. 
We begin with the IN transition. In this case, we take 



6 — aQ 



£ = (i ~ a £ ) 



(5.6) 



where Q is the value of Q at the time of crosslinking. 
We could have taken g rather than g~ x proportional to 
Q. Since we are interested in small Q, there is little 
difference between the two choices. Our goal is to recast 
/ in terms of the left strai n te nsor v to obtain a free 
energy of the form of Eq. (4.1). We begin by finding 
the equilibrium strain tensor A Q when Q = 0. Since 

there is no anisotropy when Q = 0, we have A_ = Agdij. 
A straightforward minimization of / with respect to Ao 
when Q = yields the equation of state 



nTA 2 + 5 r ( 7o -lbo = 0, 



(5.7) 



where 70 = (det A Q A^/ det g ) 1 / 2 . In the incompressible 
limit B r — > 00, this yields Aq = (det<? ) x / 2 . Setting 



=0 

AA T = Aq(S + 2v) and expanding / in powers of v using 

det 1/2 [I + 2g] = exp[iTr ln(l + 2g)} 

= 1 + Trg - [Trg 2 - ±(Trg) 2 } + • • • , (5.8) 

we obtain 

Sf = (iTrg 2 + ±B(Trg) 2 

—a/iTrvQ — ^nTTrln(S — aQ) (5-9) 

to harmonic order in v, where \i — tiT(Iq/1)A^ and B = 
B r j 2 — \i. This energy is identical to f c \ + fc of Eq. (4.1) 



plus a part depending on Q alone, which can be absorbed 



into f'q [Eq. ( f4.2[) 1. The strain can be integrated out to 
yield 



Sf 



1 



nT[(l /l)A 2 - l)]TrQ 2 + 0(Q 4 ) 



a 2 nT 



lo nT 
~T~B 



nT 

~B I "=o 



TrQ 2 . 



(5.10) 



The final form of this equat ion was obtained by solving 
the equation of state [Eq. (5.7)] for Ao to lowest order 
in nT / B and a 2 . In the incompressible limit (B = 00) 
when Iq/1 = 1, there is no shift in the coefficient of TrQ 2 
and, thus, no shift in the limit of metastability of the 
isotropic phase when the system is crosslinked in the 
isotropic phase, but there is a small shift proportional 
to TrQ 2 when it is crosslinked in the nematic phase [^6) . 

If the system is compressible, B ^ 00, or if the funda- 
mental step lengths Iq and I are different, then there is a 
shift in the coefficient of TrQ 2 even when the system is 

crosslinked in the isotropic phase. 

Deep in the nematic phase, biaxial fluctuations are 
suppressed. If we assume they are completely frozen out, 
then the step-length tensor depends only on the director, 
and we can take 
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i: 



+ (r- 1 - 1)71^} (5.11) 

and loij — lo±[Sij + (r— l)rioinoj], where r = lu/lx- Since 
^ has been scaled away by the transformation from A r 
to A, the equilibrium strain A for a given n will have 
components paral lel an d perpendicular to n and will have 
the form of Eq. ( |3.2b| ) with riQi replaced by m. As n 



rotates so does A Q , but the magnitudes A || and A j_ do 
not change. Setting 8v — and minimizing over A 
find the equations of state 

„T^A2|| +S r ( 7o - l)7o = 0, 
nT^AL +B r ( l0 - l)7o = 0. 



(5.12) 



for all B r 



These equations imply A^/A 2 ^ = r = l\\/l± 

Using det(| + + 2Sg) = det A Aq det(5 + A^lA" 1 ) 
and expanding in Sv, we obtain 



*/™=M(Hl+«) + 2 B Hl 
+2/i'ti;f| , iwfl 



isotropic phases. In both cases, the equilibrium phase 
exhibits the soft elasticity characteristic of spontaneous 
breaking of the rotational symmetry of the isotropic 
state. Thus, additional physics must be added to the 
simple neoclassical model to produce the expected mem- 
ory of the anisotropy of the nematic state at crosslinking 
and the concomitant destruction of soft elasticity. There 
are a number of mechanisms that will produce this mem- 
ory. For the purposes of illustration, we will consider 
here only a simple model studied by Verwey and Warner 
p5[ in which soft elasticity is destroyed via randomness 
in the sequence of rigid and flexible units along polymer 
chain segments. The free energy of this model reduc es as 
expected to the general form discussed in Sec. IV C . 

The sequence randomness along the chain causes the 
coupling parameter a to be a random variable with av- 
erage (a) and variance (((5a) 2 ). The chain energy [Eq. 
(|5.3| )] must be averaged over a, which appears in both g 

and g . This average (ignoring the det terms) is 



where fi = nT(l 
-l 



\\x w \\±> 
2 B 



OJ-A)||)Ao-L> 



Brit 



(5.13) 
fx, and // = 



ch/ 



InTGo/OTrA^AjV 1 } + Sf ch , (5.15) 



\n{2 + r + r- 1 ). This is identical to Eq. J3.13Q with 



where 



C l = B + 2fi, C 2 = B, C 3 = B C A = (j, and C 5 = 2// 



Sf ch = |(/ /O^Tr[(A o AV 1 ) - A (5 JAT^" 1 ) 



The free energy Sf' w of Eq. ( p,13| ) has a higher symme- 
try than the most general free energy 8f w of Eq. (4.22): 
it has only three rather than the five independent elas- 
tic constants. As a result, certain distortions will have 
the same energy in the model that do not have the same 
energy in the most general model. For example, purely 
dilational and compressional strains with A xx and A zz 
interchanged will have the same e nerg y in 8f' w but not in 
Sf w . The si mpli fied form of Eq. (5.13) resulted from our 



±nT(Z /0((M 2 )TrA QoA! 



(5.16) 



-1/2 



We can now proceed as before. Let A = A(.9 
express A m terms of v, and expand in powers in Q 
The result is 



use of Eq. ( |5.11 ) for l^ 1 . In general, l^ 1 depends on the 



(/ch) = -(l /l)nT(a)Tryg 
-±nT((fc) 2 )TrA(g )' 1/2 Q A T {g T )Q - 



(5.17) 



full tensor order parameter Qi 



S(niTij 



i<5- 

3 y 



Bi 



where B t] is the biaxial part of Q i3 with components Thi s energy is identical to Eq. ( f4.31[ ) which we expected 
in the plane perpendicular to n. Deep in the nematic 
phase, fluctuations SS in the magnitude of S and in Bij 
are small. The most general form of to lowest order 
in 8S and B^ is 



on general grounds. In the absence of strain, the second 
term of this equation tends to align the principal axis of 
Q along Q^. 



17^ + {r- 1 



l)niTij + a5S(nirij 



i%) + !>!!,,[ 



(5.14) 



where a and b are numbers. The nematic energy has 
contributions ^As(SS) 2 + ^AbTtB 2 in addition to the 
Frank free energy. Integrating out SS and B^ from the 
total free energy will yield an elastic energy in Wn, w" fc , 
and w;n_|_ with five independent elastic constants, whose 
calculation we leave to the reader. 



A. Crosslinking in the Nematic Phase 

There is no qualitative distinction in the simple neo- 
classical theory between crosslinking in the nematic and 



VI. CONCLUSIONS AND FUTURE DIRECTIONS 

In this mostly pedagogical paper we have formulated a 
classical elasticity theory of nematic liquid-crystal gels, 
carefully incorporating all underlying symmetries and 
emphasizing the distinction between independent target 
and reference space rotational symmetries. Our formu- 
lation leads to a straightforward demonstration of the 
soft elasticity of nematic-gel phases that form via spon- 
taneous symmetry breaking from an isotropic gel. This 
soft elasticity is characterized by the symmetry-enforced 
vanishing of a shear modulus and vanishing stress up to 
critical values of the appropriately applied strain. These 
and other predictions that emerge from our formulation 
are consistent with earlier predictions of the neo-classical 
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liquid-crystal rubber theory H- [0, which had been very 
successful in explaining many beautiful experiments on 
liquid-crystal elastomers. 

The advantage of our formalism is that it elucidates the 
origin of the novel soft elasticity of nematic gels, show- 
ing that it is dictated by general symmetry principles 
common to any spontaneously uniaxially ordered elas- 
tic medium and is not limited to any specific model of 
such materials. Thereby, our analysis also demonstrates 
a close connection between nematically ordered elas- 
tomers to other well-studied "soft" lattices, such as smec- 
tics (which by symmetry include cholesterics), colum- 
nar phases of fluid liquid crystals, and tensionlcss mem- 
branes, where rotational symmetry (corresponding to an 
arbitrary choice of smectic layers, columns, and mem- 
brane normal orientations) similarly enforces the vanish- 
ing of specific elastic moduli. This connection allows us 
to carry over much of the insight from those systems 
to gels. For example it seems likely that the buckling 
instability J55| in smectic liquid crystals under exten- 
sional strain parallel to layer normals will provide in- 
sight into the stripe instability |L9| of a nematic elas- 
tomer subjected to extensional strain perpendicular to 
its anisotropy axis or to the as yet unstudied generaliza- 
tion of this phenomena to compressional strain parallel 
to the anisotropy axis. 

Our formulation also permits a straightforward incor- 
poration of a variety of important effects such as spatial 
variations due, e.g., to boundary conditions, ever present 
thermal fluctuations (57]] , and local gel heterogeneity , 
thereby allowing a full statistical-mechanical treatment 
of nematic elastomers. Again, experience with smectics 
p7| , |3l| , columnar phases |H| of conventional liquid crys- 
tals, and the flat phase of tensionless elastic membranes 
pl| p>3[ strongly suggests that the latter two effects will 
qualitatively modify long scales elastic properties of ne- 
matically ordered gels, leading to phenomena such as, 
for example, anomalous elasticity, negative Poisson ra- 
tio, and topological glass order. A connection of liquid- 
crystal gels to a large body of work on closely-related 
systems of conventional liquid crystals confined in rigid 

1 1 , 3~l|] , naturally leads to an 
What role does gel elas- 



gels, such as the aerogel 
important general question 



ticity play in determining the properties and stability 
of liquid-crystal phases confined inside flexible (as op- 
posed to aerogels) heterogenous gels, such as e.g., aerosils 
jll],[l2]]? With the elastic formulation presented here we 
plan to address this question in a future publication. 

Finally, the presented description is also natural for 
treatment of fluctuating nematic elastomers membranes 
p6|] , which constitute a new universality class of mem- 
branes, adding to the well-studied classes of fluid, hexatic 
and crystalline membranes p5| . In addition to the rich- 
ness exhibited by those systems, we expect new physics 
associated with the interplay of the in-plane and undu- 
lation nonlinear elasticity, both expected to be impor- 
tant in elastomer membranes |26[ ] . Finally, such in-plane 
orientationally-ordered elastic membranes are novel real- 



izations of anisotropic membranes, predicted to exhibit 
flat, tubule, and crumpled phases paL subsequently ob- 
served in Monte Carlo simulations |59|]. We plan to ex- 
plore these and other phenomena and realizations of ori- 
cntationally ordered elastomers in future publications. 

Pel 
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APPENDIX A: REVIEW OF NEMATIC ENERGY 

In this appendix, we will review standard treatments of 
the isotropic-to-nematic transition, principally to estab- 
lish notation. We introduce a complete set of orthonor- 
mal symmetric-traceless tensors satisfying Trl_ a lf = 





4 o 

1 



1 



1 
_ . -1 
v2 \ o 



I 2 = 



J 4 



1 



1 
_ . 1 
v2 \ o 




1 

71 




1 
1 



(Al) 



Any symmetric-traceless tensor can be expressed as a 
linear combination of these matrices: Qij = ^2 a=0 Qal"j, 
where Q a = TrQI_ a . Thus, 



Qo — y z[Qzz 
Q2 = V2Q xyi 



Q 3 = V2Q XZ , 



Qi 



V2Q yz . 



°iyyii 

(A2) 



/372>°n° - (1/3)^). 



With n° along the z axis, if, — y «j/ ^ V "i '"j — 
In uniaxial nematic phase, = 5*(n°n° — (l/3)<5y) = 

Qolfj and Qo = \/2/35. The Landau-de Gennes free 
energy for a nematic is 



fa = ±r Q Tr£ 2 - w 3 Tvg + w 4 (TrQ 2 ) 2 
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-w 3 



2 V 3^0 



QoiQl + Ql) 



^w3[ZQi{Ql - Ql) + 6Q2Q3Q4 



2^0(01 + Qt) 
(A3) 



Minimization with respect to Qo yields the equation of 
state 



r Q Q + 4w 4 Qo - V ! w 3<2o = 



(A4) 



— 4/2 
B 1 =5 1 + - Mo -2^- 7 

— 2 1 [2 



_ 1 

-B 3 = Si - /i 2 + -/x 

Ba = /i 2 . 



Ti' 



(B4) 



Then expansion to second order in deviations 5Q a 
Q a ~ Q° a yields 



5f=U 1 {8Sf+ 1 MQ\ + Ql) 



(A5) 



APPENDIX C: LINEARIZED LIMITS OF 
EULERIAN AND LAGRANGIAN ELASTICITY 



where 



M = l{r Q + 8 Wi S 2 -2w 3 S) 
A 2 = 3w 3 S 



(A6) 



As anticipated from the underlying rotational invariance, 



there are no terms proportional to 5Q\ or 5Q\ 



APPENDIX B: EVALUATION OF B A 

In this a ppend ix, w e out line the algebraic steps be- 
tween Eqs. ( 4.11 ) and ( 4.12 ). We need to integrate over 
SS and Q\ and Q 2 . Since these variables appear only to 
quadratic order, the integration is trivial and yields 



6f u _ v = iBi(TMu) 2 +p 2 (t;?- 
+HqSvq — jSvqTtSu, 



(Bl) 



It is often the case that a linearized theory of elastic- 
ity, in which nonlinear strains are replaced by their lin- 
earized limits and only terms to harmonic order in these 
linearized strains are included in the free energy, pro- 
vides an adequate description of elastic distortions. It is, 
therefore, interesting to see how this linearized limit is 
reached. It turns out that this limit can be taken more 
cleanly in the Eulerian picture in which the displacement 
field is a part of the phase of a mass-density wave rather 
than the Lagrangian picture in which u(x) is a displace- 
ment relative to a reference configuration. For a further 
discussion of these two pictures of elasticity, see Ref. [pi . 
Much of our intuition about how to construct a linearize 
theory comes from the Eulerian picture in which the dis- 
placement field is a vector field in space that obeys the 
usual rules of transformation of vector fields. In this Ap- 
pendix, we will discuss the linearized limits of Eulerian 
and Lagrangian elasticities. 



where 



Bi = B 



Mo = M " 



^2 = 



7 = 



where 



16s 2 S 2 

9A' 
4^ 



a 2 + (2tyfi) 



- 4, 3 



2\ 3/2 stS 



A' 



A'=A 1 + %S 2 + 3 —. 



(B2) 



(B3) 



Setting Trim = (5v xx + Sv yy + Sv zz ), replacing v , v\, 



and V2 with expressions obtained from Eqs. (A2) (with 
the tensor Q replaced by Sv) and including the Sv a — 



(t/fj,) SQ a with a = 3,4 terms in Eq. (4.11), we obtain 
Eq. (§T2j) with 



1. Eulerian Elasticity 

In Eulerian elasticity, the displacement field u(x) is 
a vector field in three space. Like all vector fields that 
transform under the same group as space itself, u trans- 
forms under a rotation of the whole sample as 



u'(x) = C/utx)^ 1 = Ou(0 _1 x) 



(CI) 



where U is a rotation operator (e.g., quantum mechanical 
operator) and O is its associated 3d rotation matrix. Here 
the prime indicates the value of the field after the rotation 
operator is applied. To leave the system unchanged, U 
must be an operation in the point group of the crystal. In 
the Eulerian picture, u(x) is a Goldstone field associated 
with the broken spatial symmetry of a crystal. Thus, 
strictly speaking the highest-symmetry point group in 
three dimensions in the cubic group. To make contact 
with our discussion of gels, we can, however imagine a 
system in which all rotations are in the point group. Since 
u is a vector field, diUj is a tensor field that satisfies 



UdiUjU 



O ik d' kUl {x!)Oi 3 \ (C2) 
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where x' = 0~ 1 x, d[ — djdx\, and as before the prime 
indicates the value of the operator after rotation. Alter- 
natively, we can introduce 

riij = djUi, (C3) 

which 

rf =gv QT 1 (C4) 

Scalars created from r] or from u and its derivatives are 
invariant under U. For example, 

Vii = diUi, T) i:j T)ji (C5) 

etc. are scalars under U. 

The above symmetries and considerations apply to any 
vector field. The displacement field u, however, has ad- 
ditional properties arising from the fact that it is a Gold- 
stone field. In particular the systems is invariant under 
rotation of the mass-density wave crystal (which is not 
the same thing as rotating the whole sample). The trans- 
formation 

x - u(x) -» £(x - u(x)) (C6) 

rotates the crystal. Thus, the transformation 

u -» u'(x) = -0x - u(x)) + x (C7) 

does not change the energy of the system. This implies 
that the elastic energy will depend only on the Eulerian 
symmetrized strain: 

ufj = diUj + diUj - d^kdjUk (C8) 
« d^j + d^j (C9) 

where the final form is its harmonic limit. 

Thus, we have two symmetries: (1) symmetries associ- 
ated with rotation of the whole sample, and (2) rotations 
of the lattice. Invariance with respect to the first requires 
that the energy depend only on scalars formed by con- 
tracting indices of both gradients and u's. The second 
invariance requires that the energy be a function only 
of the nonlinear strain ufj . The interesting thing is that 
both uf- and its linearized form transform like tensors un- 
der (1), i.e., under O. Thus, contracted tensors of either 
u E or its linearized form are scalars under O. 



Under infinitesimal rotations, 

Onj = 5 i3 + ei ok O T k (Cll) 
Rij = S i:j + e iok 9 Rk . (C12) 

and 

diu'j = fyiij + e ijp (6 Rp - 6 Tp ) 

+tikpd' k u j 9 Rp - ekjpd'jUkOTp (C13) 

The energy is invariant under independent rotations 
through 9 R and 9t- As we have seen, these invariances 
are guaranteed by making the free energy a function only 
of the fully contracted ity or Vij tensors. 

Now, let us look at the linearized limit. Under both 
rotations, we have, 

diu'j + dju'i = d[uj + djUi 

+ (tikpd' k Uj + ejkpd' k Ui)6 Rp 
-{ekjp&iUk + e k i V d' u k )Q T p- (C14) 

Note that to leading order in u and 9, this symmetrized 
combination is independent of 9 R and 9t as it should be. 
The terms of order 9u tell us about the tensorial rotation 
properties of the system. // Q R = Q T = 9, then the 
symmetrized combination uf^ = diUj + djUi transforms 
like a tensor, i.e., 

^^Q^Qr 1 . (C15) 

The linearized strain u S does not, however, transform 
like a tensor under independent rotations 9 R and 9t- If, 
for example, 9 R is zero, the term proportional to the 
product 9 Ri u k in the transformation of u s depends on 
both the symmetric and anti-symmetric parts of diUj. 
Only the fully nonlinear strains tiy and transform like 
tensors even to linear order in 9 R or 9t- We leave it as 
an exercise to verify this explicitly. Thus, we cannot use 
the linearized tensors to discuss the rotation and tensorial 
properties of strains in the Lagrangian language. We can, 
however, as discussed above use them fruitfully in the 
Eulerian language. However, once we have constructed a 
rotationally invariant Lagrangian energy, we can replace 
nonlinear strains by linearized ones to discuss harmonic 
clastic fluctuations. 



2. Lagrangian Elasticity 

In Lagrangian Elasticity, there are as we have discussed 
two symmetries: (1) rotations Or in the target space 
and (2) rotation O r in the reference space. Under these 
operations, the displacement vector satisfies 

R'(x) =fi T R(fi fl 1 x) (CIO) 
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